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NUMERICAL  SIMULATIONS  OF  THE  CIRCULATION  OF  THE  EASTERN  MEDITERRANEAN 


by 

Despina  Grillaki  and  Steve  Piacsek 


ABSTRACT 


Oceanographic  observations  have  shown  that  the  surface  circulation  in  the 
eastern  Mediterranean  consists  of  the  zonal  North  African  Current  along 
with  systems  of  gyres  and  rings  both  north  and  south  of  the  Current.  A 
large  cyclonic  gyre  occurs  between  Cyprus  and  Crete,  and  another  between 
Greece  and  Libya.  A  weaker,  quasi -permanent  eddy  is  formed  in  the  Ionian 
sea  and  is  intensified  by  the  flow  through  the  Strait  of  Sicily.  An 
anticyclonic  ring  occurs  south  of  the  North  African  Current  in  the  Gulf  of 
Sirte  and  another  occurs  in  the  southwest  Levantine  Sea.  The  reduced- 
gravity  model  (one  active  layer)  of  Hurlburt  and  Thompson  <6>  has  been 
adapted  and  applied  to  the  semi -enclosed  basin  of  the  eastern 
Mediterranean.  Using  monthly  averages  of  climatological  wind  stresses  and 
the  estimate  of  the  annual  mean  flow  through  the  Strait  of  Sicily, 
numerical  simulations  of  the  circulation  in  the  eastern  Mediterranean  basin 
have  been  carried  out.  Model  runs  both  with  and  without  inflow  through  the 
Straits  have  been  performed.  The  overall  results  agree  with  the 

observations.  They  show  the  presence  of  the  North  African  Current,  two 

persistent  cyclonic  gyres  of  varying  strengths  and  two  anticyclonic  rings. 

This  modelling  effort  indicates  the  importance  of  the  wind  induced 
circulation  versus  transport  for  the  east-central  part  of  the  eastern 
Mediterranean.  It  is  a  first  step  towards  a  combined  hydrodynamic- 

thermodynamic  model  of  the  circulation  in  that  area. 
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INTRODUCTION 

> 

This  SACLANTCEN  report  describes  a  modelling  effort  that  uses  realistic 
'  geometry,  combined  forcings,  mass  transport,  and  climatology  to  shed  some 

[  light  on  the  circulation  of  the  eastern  Mediterranean  and  its  driving 

,  mechanisms.  In  particular  it  discusses  the  relative  importance  of  wind 

i  forcing  vis-a-vis  inflow-forcing  through  the  Strait  of  Sicily.  The 

'  hydrodynamic  numerical  model  simulates  the  major  observed  hydrodynamic 

features  of  the  upper  layer  circulation,  such  as  meandering  currents, 
persisting  gyres  and  transient  eddies. 

•  The  Strait  of  Sicily  separates  the  Medi terranean  into  western  and  eastern 

,  basins  (Fig.  la).  The  eastern  Mediterranean  basin  is  subdivided  into  the 

Adriatic,  the  Ionian,  the  Levantine  and  the  Aegean  Seas. 

The  surface  circulation  in  the  eastern  Mediterranean  (Fig.  lb)  is 
significantly  influenced  by  both  the  Atlantic  inflow  (originally  entering 
the  Mediterranean  through  the  Strait  of  Gibraltar)  and  the  Levantine  basin 
water.  The  Atlantic  water  enters  the  eastern  Mediterranean  through  the 
Strait  of  Sicily  (142  km  wide,  400  m  deep)  and  forms  a  slow  eastward- 

|  flowing  surface  layer,  overriding  the  sinking,  more  saline  Levantine  water 

that  moves  westward  as  an  intermediate  lower  layer  <8>. 

The  Atlantic  water  has  been  traced  into  the  eastern  Mediterranean  by  summer 

observations  of  a  salinity  minimum  in  the  upper  25  to  75  m  <8>.  Modified 

progressively  by  evaporation  of  the  surface,  the  Atlantic  water  transits 
the  Strait  of  Sicily,  flows  along  the  Libyan  and  Egyptian  coasts  and  turns 
along  the  Israeli  coast  into  the  northeast  Levantine  sea  (south  of  Crete). 
A  series  of  cyclonic  eddies  form  north  of  the  Current. 

A  smaller  branch  of  the  Atlantic  inflow  through  the  Strait  of  Sicily 

1  deflects  to  the  North  and  is  traced  into  the  Ionian  sea. 

The  prevailing  westerly  winds  favour  the  above  circulation  pattern.  These 
perennial,  relatively  strong  winds,  acting  on  the  limited  area  of  the 
Mediterranean  basin,  generate  large  eddies.  These  eddies  are  associated 
with  vertical  motions  and  play  an  important  role  in  water  mass  formation  by 
bringing  underlying  water  parcels  closer  to  the  mixing  regions  <4>. 

Meteorological  forcing  and  thermodynamics  are  the  main  driving  physical 
mechanisms  of  the  eastern  Mediterranean.  With  the  help  of  a  simple 
hydrodynamic  model,  we  have  studied  the  circulation  and  dynamics  of  the 
eastern  Mediterranean. 

The  model  consisted  of  2-layers,  an  upper  layer  200  m  deep  and  a  lower 
layer  of  infinite  depth.  This  'reduced  gravity'  model  has  been 
successfully  used  in  recent  years  in  several  studios  about  the  circulation 
and  dynamics  of  semi -enclosed  seas  such  as  the  Gulf  of  Mexico. 
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The  recent  advances  in  computer  technology  have  marie  it  possible  to  run 
numerical  ocean  models  with  sufficient  resolution  to  depict  mesoscale  ocean 
features.  The  most  extensive  study  along  these  lines  was  done  by  Hurlburt 
and  Thompson  <6>,  in  which  they  successfully  simulated  many  dynamic 
features  of  the  loop  current  in  the  Gulf  of  Mexico. 

More  relevant  to  this  work  in  the  Mediterranean  were  the  modelling  studies 
of  Preller  and  Hurlburt  <  1 2 >  and  Preller  <13>  for  the  Alboran  Sea  and 
Preller  and  Heburn  <14>  for  the  western  Mediterranean  Sea.  The  western 
Medi terranean  study  also  used  realistic  geometry  and  bottom  topography,  and 
combined  forcings  from  the  Strait  of  Gibraltar  inflow  and  monthly  averaged 
wind  stresses. 

As  for  the  eastern  Mediterranean,  one  earlier  study  <]4>  used  the  reduced- 
gravity  model  with  an  idealized  rectangular  basin  geometry.  Another 
eastern  Mediterranean  study  (unpublished)  uses  a  two-layer  active  'rigid 
lid'  model  (i.e.  assumes  a  level  surface).  Both  of  these  studies 
simulated  only  the  wind-driven  circulation. 


FIG.  1  a)  THE  MEDITERRANEAN  AND  ITS  MAJOR  SEAS. 

b)  SURFACE  CIRCULATION  OF  THE  MEDITERRANEAN , 
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1  THE  MODEL 


The  numerical  model  used  in  this  study  of  the  eastern  Mediterranean  has 
been  modified  and  adapted  from  one  developed  for  the  Gulf  of  Mexico  <6>. 

The  model  is  an  active  upper  layer,  hydrodynamic  model  on  a  b-plane.  The 
model,  of  the  first  baroclinic  mode,  exclude  baroclinic  instability  and 
assumes  Boussinesq  approximation.  Such  models  are  known  as  'reduced- 
gravity  models'  because  the  effective  gravity  acting  on  the  active  upper 
layer  is  reduced  by  the  density  ratio. 


g'  =  g(Ap  j  p  )  (Eq.  1) 

and 

Pl  -  P2  (  p  .  (Eq.  2) 


where 

g  is  acceleration  due  to  gravity  and 
p  is  the  density. 


The  active  upper  layer  is  200  m  deep  and  the  static  lower  layer  is 
infinitely  deep;  between  them  is  an  impermeable  interface.  The  mode 1  uses 
realistic  irregular  geometry  (coastline),  Fig.  2,  drawn  from  the  Synthetic 
Bathymetric  Profiling  System  (SYNBAPS)  of  Vanwyckhouse  <15,  16>. 


DX.DY  =  0.250,0.167  (DEG)  DBT  =1000.0  [Ml 


FIG .  2  BATHYMETRY  ON  THE  EASTERN  MEDITERRANEAN  SEA 
(from  SYNBARS  data). 
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The  grid  resolution  is  1/4  of  a  degree  (21.4  km  in  the  east-west  direction) 
and  1/6  of  a  degree  (18.6  km  in  the  north-south  direction).  The  fine  grid 
resolution  allows  the  study  of  mesoscale  features. 

The  model  domain  extends  from  10°  to  36°E  longitude  and  from  30°  to  38°N 
latitude.  For  the  coastline  a  cut-off  depth  of  20  m  was  used  in  order  to 
eliminate  errors  existing  at  the  boundary  of  the  data  set.  Except  at  the 
ports,  the  boundaries  are  rigid  with  a  no-slip  boundary  condition.  The 
inflow  transport  values  are  based  on  observations. 

The  vertically  integrated  equations  used  in  the  model  are  the  equation  of 
motion  (Eq.  1),  including  pressure,  Coriolis  acceleration,  wind  stress  and 
frictional  terms,  and  the  equation  of  continuity  (Eq.  2).  They  are: 


eVi 


—  +  (V'Vi  +  Vx  -  V)  vi  +  K  x  f  Vi 


-  hi  !  P|  +  p  (T1-T2)  +  A 

(Eq.  3) 


ehi 

et 


+  v*v.  =0 


(Eq.  4) 


where 


v  =  —  i  +  ®_  j  , 

0  x  0y 


pl  =  9  ►>!  , 


vi  =  h[  vi  =  hi  (uj  i  +  V[  j) 


91  =  9(p2 -Pi )  p' 


and 


f  =  ^0  +  B(y-y0)  , 


Ti  =  Ti  (x)  i  +  Ti  (y)  j 


x  and  y  are  tangent-plane  Cartesian  coordinates,  U[  and  vi  are  the 
eastward  and  northward  velocity  components  in  the  upper  layer,  ht  is  the 
upper  layer  thickness,  t  is  the  time,  g  is  the  acceleration  due  to 
gravity,  pi  is  the  density  of  seawater  in  layer  i  ,  f0  and  V0  are 
the  values  of  the  Coriolis  force,  and  the  y-coordinate  at  the  southern 
boundary,  Tx  is  the  wind  stress,  and  Ty  is  the  interfacial  stress.  The 
remaining  parameters  are  defined  in  Table  1. 
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TABLE  1 

MODEL  PARAMETERS 


PARAMETER 

DEFINITION 

VALUE 

A 

Eddy  viscosity 

2  -1 

100D  m  s 

b 

(Df/Dy) 

.11  1  -1 

1.8  x  10  ms 

fo 

Coriolis  parameter 

5  .1 

8  x  10  s 

9' 

Reduced  gravity  due  to 
stratification 

0.02  m  s'2 

hi 

Undisturbed  upper  layer 
depth 

200  m 

Lx 

East-west  dimension  of 
the  model  domain 

10°E  to  36°E  longitude 

Ly 

North-south  dimension 
of  the  model  domain 

30°N  to  38°N  latitude 

x  .  y 

Horizontal  grid  reso¬ 
lution 

21.4  km  (1/4°)  x 

18.6  km  (1/6°) 

t 

Time  step 

2700  s 

Vin 

Inflow  velocity 

5  cm  s'' 

Tsin 

Inflow  spin  up  time 
constant 

30  days 

^swind 

Wind  spin  up  time 
constant 

2  days 

a 

Angle  of  inf’ow 
(Strait  of  Sicily) 

1 35  0 

Inflow  is  exactly  balanced  by  outflow  through  an  open  eastern  boundary 
(port).  This  is  accomplished  by  allowing  normal  flow  at  the  eastern 
boundary  to  be  self-determined  and  by  imposing  an  integral  constraint  on 
total  mass  outflow. 

The  model  parameters  and  their  values  are  given  in  Table  1.  The  value  for 
the  eddy  viscosity  was  chosen  to  keep  the  Reynolds  number  low  for 
computational  stability.  The  value  of  the  reduced  gravity  coefficient  is 
the  average  density  difference  of  the  Mediterranean  environment  at  these 
levels.  The  undisturbed  upper  layer  depth  value  generally  represents  the 
upper  end  of  the  thermocline  and  was  derived  from  observations. 
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The  inflow  velocity  is  spun  up  over  a  period  of  thirty  days  to  minimize  the 
excitation  of  high  frequency  waves.  The  winds  have  a  shorter  spin  up  period 
of  two  days. 

The  inflow  port  representing  the  Strait  of  Sicily  is  defined  by  5  grid 
points  and  is  142  km  wide.  The  artificial  outflow  port  accounting  for  the 
evaporation  and  sinking  of  waters  into  greater  depths  is  also  defined  by 
S  grid  points  but  is  93  km  wide.  An  inflow  velocity  of  5  cm/s  corresponds 

to  a  transport  of  0.96  x  106  m3 /s  (0.96  Sv)  through  the  Straits  of  Sicily, 
which  agrees  with  observations  <2,  3>. 

The  direction  of  the  inflow  is  taken  as  135°.  The  inflow  is  independent  of 
time  and  represents  only  the  mean  annual  flux.  Its  temporal  and  spatial 
variations  in  the  Strait  of  Sicily  are  assumed  to  balance  out,  i.e.  no  net 
change.  However,  in  reality,  the  Atlantic  water  travels  by  spreading 
laterally  and  is  seasonally  dependent. 

The  open  eastern  boundary  condition  outflow  port  accounts  for  the 
evaporation  and  the  sinking  of  surface  waters  and  therefore  satisfies  the 
equation  of  continuity. 

The  outflow  port  location  was  chosen  near  the  Israeli  coast  as  a  first 
attempt  to  approach  the  observed  oceanographic  circulation  patterns. The 

model  results  have  shown  that  the  outflow  port  location  was  not  an 
important  factor  for  this  study.  The  Atlantic  origin  water  sinks  in  the 

northeast  Levantine  Sea  and,  in  the  Crete-Rhodes-Cyprus  area,  is 
transformed  into  the  denser  Intermediate  Medi terranean  water.  In  winter 

the  surface  atmospheric  fluxes  enhance  the  vertical  circulation  which 
accompanies  the  formation  of  the  Intermediate  Mediterranean  water. 

The  monthly  climatological  wind  forcing  has  been  derived  from  twenty  years 
of  ship  observations  <9>.  These  monthly  averages  of  wind  stresses  are 

computed  on  a  1°  latitude  x  1°  longitude  (Marsden  Square)  grid.  When 
applied  to  the  1/4°  by  1/6°  model  grid,  the  winds  are  bilinearly 
interpolated. 
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2  MODEL  RESULTS 


The  model  was  subjected  to  a  series  of  experiments  to  determine  the 
influence  of  mass  transport  and  wind  forcings,  singly  or  in  combination,  on 
surface  circulation. 

Firstly,  the  model  was  forced  only  by  mass  transport  through  the  Strait  of 
Sicily  to  examine  the  contribution  to  the  surface  circulation.  Secondly, 
the  model  was  forced  by  monthly  climatological  wind  stress  to  examine  the 
purely  wind-driven  circulation.  Thirdly,  a  combination  of  transport  and 
wind  stress  forcings  were  applied. 


2.1  Mass  transport  (Inflow)  forcing  experiment 


The  experiment  has  been  integrated  over  five  years.  The  results  are 
presented  (Figs.  3,  4)  in  terms  of  pycnocline  height  deviation  of  the 
interface  between  the  upper  and  lower  layers  and  current  patterns.  (The 
pycnocline  height  deviation  is  the  deviation  of  the  interface  from  its  rest 
position  of  200  m).  Figure  3  presents  results  for  January,  February,  March 
and  April  for  the  5*^  year  of  simulation.  Figure  4  shows  the  current 
pattern  of  the  active  upper  layer  for  the  month  of  December, 

With  inflow  forcing  alone,  the  Atlantic  origin  water  passes  through  the 
Strait  of  Sicily,  and  flows  along  the  African  coast  as  an  extension  of  the 
North  African  Current.  Little  difference  is  observed  between  figures,  thus 
indicating  the  existence  of  a  steady  state. 

No  seasonal  changes  are  found  in  the  circulation  patterns.  However,  the 
effects  of  changes  in  direction  and  velocity  of  the  inflow  do  affect  the 
path  of  the  current,  as  in  the  Alboran  sea  <12>  but  this  effect  was  not 
studied.  The  flow-channel  of  the  Strait  of  Sicily  is  oriented  such  that 
the  entering  current  flows  southward  directing  the  mass  transport  close  to 
the  North  African  coast. 

Compared  to  observational  data  <10,  1 1  > ,  the  model  indicates  only  the 
North  African  Current  with  very  low  speeds  and  an  absence  of  gyres  and  of 
circulation  in  the  northern  part  of  the  basin  (Fig.  5).  Because  the  inflow 
alone  represents  only  a  minor  contribution  to  the  east  Mediterranean 
circulation,  the  variations  in  inflow  and  the  outflow  port  location  were 
not  important  for  this  particular  study. 


2.2  Wind  forcing  experiment 

In  this  experiment  the  ports  are  closed  and  the  wind  stresses  (analyzed  by 
May  <9>)  blow  over  the  eastern  part  of  the  basin.  Thus  the  solutions 
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(Fig.  6)  represent  circulation  due  only  to  climatological  wind  stress 
applied  in  a  closed  basin.  Figures  7,  8,  9  show  monthly  solutions  for  the 

fifth  year  of  integration  in  terms  of  the  pycnocline  anomaly  and  Fig.  10 

shows  the  current  pattern  of  the  active  upper  layer  for  the  month  of 

December. 

In  the  central  part  of  the  eastern  basin  several  persistent  gyres  appear 
along  with  a  strong  cyclonic  circulation.  Cyclonic  gyres  persist  between 
southern  Greece  and  Libya  and  between  Cyprus  and  Crete.  Also  a  permanent 
anticyclonic  gyre  appears  in  the  southwest  Levantine  basin. 

For  the  western  part  of  the  eastern  Mediterranean  basin  the  circulation 
pattern  varies  seasonally.  The  North  African  Current  appears  in  the 
solution.  A  cyclonic  quasi -steady  eddy  appears  in  the  Ionian  Sea  and  a 

seasonal  anticyclonic  gyre  appears  in  the  Gulf  of  Sirte. 

Model  solutions  show  monthly  and  seasonal  variation  of  circulation  with  an 
maximum  in  winter  and  a  minimum  in  late  summer.  However,  the  Aegean  is  an 
exception;  here  the  maximum  occurs  in  late  summer.  The  magnitudes  of  the 
wind  stresses  (Fig.  6)  also  show  an  overall  maximum  in  the  winter  with  the 
exception  of  the  Aegean  Sea,  where  the  maximum  occurs  in  late  summer. 
(These  strong  northern  summer  winds  are  called  the  Etesians.)  May  <9>  has 
determined  large  wind  stress  values  in  the  area  of  the  two  major  cyclonic 
gyres  formed  east  and  west  of  Crete. 

The  model  results  compare  well  with  the  observational  data  (Fig.  5)  in  the 
central  and  eastern  parts  of  the  basin  but  give  less  consistent  all  year 
around  results  for  the  western  part.  The  results  indicate  that  the  wind 
stress,  as  a  driving  force,  is  responsible  for  the  development  of  the  major 
cyclonic  gyres  throughout  the  year. 


2. 3  Combined  forcing  experiment 

This  experiment  consisted  of  combining  the  mass  transport  through  the 
Strait  of  Sicily  with  the  wind  driven  experiment  at  the  end  of  the  fifth 
year.  Then  this  combination  was  run  for  five  more  years.  This  approach 
was  used  because  the  inflow-forced  solutions  reach  a  steady  state  more 
quickly  than  the  wind-forced  ones. 

Figures  11  to  14  represent  the  model  solution  and  resemble  a  superposition 
of  the  two  previously  examined  independent  cases.  Several  results  are 
apparent:  the  inflow  seems  to  intensify  both  the  semi -permanent  eddy  in 
the  northwestern  basin  and  the  anticyclonic  ring  in  the  Gulf  of  Sirte, 
which  now  becomes  permanent  throughout  the  year.  The  North  African  Current 
appears  at  the  south  of  the  basin.  In  the  eastern  basin  the  cyclonic  gyres 
between  Cyprus  and  Crete  and  between  Greece  and  Libya  are  still  maintained 
as  is  the  anticyclonic  ring  in  the  southwest  Levantine  Sea.  The  intensity 
of  the  features  seems  to  grow  stronger  in  late  winter  and  early  spring. 

The  model  agrees  well  with  Ovchinnikov's  observations  <10,  11>.  Both  the 
observations  and  the  model  present  cyclonic  gyres  north  of  the 
North  African  Current  and  anticyclonic  rings  south  of  it. 
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The  similarity  between  the  modelled  (see  Fig.  5)  and  the  observed 
circulation  patterns  is  stronger  for  the  eastern  and  central  part  of  the 
basin  than  for  the  western  part.  The  implication  is  that  the  wind,  rather 
than  inflow,  is  the  dominant  physical  influence  on  the  surface  patterns. 

For  the  western  part  of  the  basin  the  restricted  domain  of  the  model  and 
the  driving  forces,  i.e.  mass  transport  and  wind  stress,  have  only  a 
limited  capability  to  describe  the  surface  circulation.  The  modelling 
experiments  also  show  that  the  inflow  represents  only  a  minor  contribution 
to  the  eastern  Mediterranean  circulation. 


SUMMARY 


The  surface  circulation  of  the  Mediterranean  east  of  the  Strait  of  Sicily 
has  been  described  by  a  simple,  one  (active)  layer,  reduced-gravity  model. 
Although  this  model  doesn't  consider  bottom  topography,  thermodynamics,  or 
mixing  (except  for  lateral  diffusion  of  momentum),  it  has  produced  a  first 
approximation  to  the  mass  transport  through  the  Strait  of  Sicily  and  the 
seasonally  varying  climatological  wind  stresses.  The  modelling  experiments 
first  determined  the  influence  of  mass  transport  on  surface  circulation. 
This  was  followed  by  experiments  to  determine  the  influence  of  wind  stress. 
Finally,  the  modelling  experiments  considered  a  combination  of  both 
factors. 


The  overall  results  showed  a  qualitative  agreement  with  Ovchinnikov 
<10,  11>  observed  winter  circulation  patterns  in  this  area.  The 

North  African  Current,  two  major  cyclonic  gyres  and  a  seasonal  eddy  to  the 
north  and  two  anticyclonic  rings  to  the  south,  observed  in  the  circulation 
patterns,  appear  also  in  the  model.  However,  they  differ  in  the  intensity 
of  circulation  and  the  location  of  the  gyre  centers.  In  the  model,  the 
large  gyres,  east  and  west  of  Crete,  are  centered  much  closer  to  the  island 
than  indicated  by  the  observed  results.  The  ring  in  the  southwest 

Levantine  basin  is  more  distant  from  the  coast  in  the  model  than  in  the 

observed  results. 

The  model  creates  a  realistic  large  scale  surface  circulation  pattern. 
Wind-induced  circulation  is  an  important  physical  mechanism  for  influence 
of  the  circulation  in  the  central  and  eastern  parts  of  the  eastern 
Mediterranean.  Mass  transport  is  only  a  minor  influence  except  in  the 

function  of  the  North  African  Current. 

This  report  presents  preliminary  results  of  experiments  to  model  the 
circulation  of  the  eastern  Mediterranean.  Additional  hydrodynamic 
computations  are  needed  to  improve  the  model.  These  should  include  the 
time  dependence  of  the  inflow  through  the  Strait  of  Sicily  and  should 

employ  a  two  (active)  layer  model  to  account  for  the  subsurface  outflow  and 
bottom  topography.  A  more  complete  numerical  simulation  of  the  circulation 
in  the  eastern  Medi tterranean  should  consider  thermodynamics  and  use  more 
realistic  wind  forcing  functions. 
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INTERFACE  DEVIATIONS 


(DEG i  36E 


a 

DAY  =  1470 
DH  =  5.0  (m) 
MIN  =  -1.05E+01 
MAX  =  1.52E+01 


b 

DAY  =  1500 
DH  =  5.0  (m) 
MIN  =  -1.06E+01 
MAX  =  1.53E+01 


38N 


i  DEC) 


3  Or  I 


C 


DAY  =  1530 
DH  =  5.0  (m) 

MIN  =  -1.04E+0I 
MAX  =  1.53E+01 


d 

DAY  =  1560 
DH  =  5.0  (m) 

MIN  =  -1.07E+01 
MAX  =  1.53E+01 


FIG,  3  SOLUTIONS  OF  THE  INFLOW  FORCES  EXPERIMENT  DURING 
THE  FIFTH  YEAR, 

(a)  January  (c)  March 

(b)  February  (d)  April 

Solid  contours  indicate  deepening  of  the  interface 
and  anticyclonic  circulation /  dashed  contours  indicate 
shallowing  and  cyclonic  circulation. 
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INTERFACE  DEVIATIONS 


a 

DAY  =  1470 
DH  =  20.0  (m) 
MIN  =  -1.31E+02 
MAX  *  4.73E+01 


b 


DAY  =  1500 
OH  =  20.0  (m) 
MIN  =  -1.20E+02 
MAX  =  4.80E+01 


C 

DAY  =  1530 
OH  =  20.0  (m) 
MIN  =  -1.16E+02 
MAX  =  4.94E+01 


d 

DAY  =  1560 
DH  =  20.0  (m) 
MIN  =  -1.16E+02 
MAX  =  4.97E+01 


SOLUTIONS  OF  WIND  FORCING  EXPERIMENT  DURING 

THE  FIFTH  YEAR. 

fa)  January  fc)  March 

(b)  February  fd)  April 

Solid  contours  indicate  deepening  of  the  interface 
and  anticyclonlc  circulation;  dashed  contours  indicate 
shallowing  and  cyclonic  circulation . 


-1.10E+02 

5.50E+01 


MIN  =  -1.11E+02 
~L  MAX  =  5.58E+01 


FIG.  8  SOLUTIONS  OF  WIND  FORCING  EXPERIMENT  DURING 
THE  FIFTH  YEAR. 

(a)  May  (c)  July 

(b)  June  (d)  August 

Solid  contours  Indicate  deepening  of  the  interface 
and  anticyclonic  circulation;  dashed  contours  indicate 
shallowing  and  cyclonic  circulation. 
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INTERFACE  DEVIATIONS 


a 


DAY  =  1710 
DH  =  20.0  (m) 
MIN  =  -1.12E+02 
MAX  =  5.46E+01 


b 


DAY  =  1740 
DH  =  20.0  (m) 
MIN  =  -1.12E+02 
MAX  =  5.28E+01 


C 


DAY  =  1770 
DH  =  20.0  (m) 
MIN  =  -1.10E+02 
MAX  =  5.32E+01 


d 


DAY  =  1800 
DH  =  20.0  (m) 
MIN  =  -1.10E+02 
MAX  =  5.09E+01 


FIG.  9  SOLUTIONS  OF  WIND  FORCING  EXPERIMENT  DURING 
THE  FIFTH  YEAR. 

(a)  September  (c)  November 

(b)  October  (d)  December 

Solid  contours  indicate  deepening  of  the  interface 
and  anticgclonic  circulation ;  dashed  contours  indicate 
shallowing  and  cyclonic  circulation. 
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oor  =  1800  usjys 


MAX  PL0TTED  SPEED  =  0.32  (M/SEC) 


FIG.  10  VELOCITY  FIELD  FOR  THE  WIND  DRIVEN  CASE. 
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INTERFACE  DEVIATIONS 


FIG.  11  SOLUTIONS  OF  COMBINED  WIND  AND  INFLOW  FORCING  EXPERIMENT 
DURING  THE  FIFTH  YEAR. 

(a)  January  (c)  March 

(b)  February  (d)  April 

Solid  contours  indicate  deepening  of  the  interface 
and  ant icyclonic  circulation/  dashed  contours  indicate 
shallowing  and  cyclonic  circulation. 


3270 

20.0  (m) 
-1  .24E+02 
6.32E+0.1 


3300 

20.0  (m) 

-1.32E+02 

6.18E+01 


3330 

20.0  (m) 
-1.35E+02 
6. 12E+01 


3360 

20.0  (m) 

-1.35E+02 

6.08E+01 


19 


SACLANTCEN  SR-92 


INTERFACE  DEVIATIONS 


a 


DAY  =  3390 
DH  =  20.0  (m) 
MIN  =  -1.31E+02 
MAX  =  6.02E+01 


b 


DAY  =  34  20 
OH  =  20.0  (m) 
MIN  =  -1.26E+02 
MAX  =  6.00E+01 


C 


DAY  =  3450 
OH  =  20.0  (m) 
MIN  =  -1.21E+02 
MAX  =  6.09E+01 


d 


DAY  =  3480 
DH  =  20.0  (ra) 
MIN  =  -1.18E+02 
MAX  =  6.24E+01 


FIG.  12  SOLUTIONS  OF  COMBINED  WIND  AND  INFLOW  FORCING  EXPERIMENT 
DURING  THE  FIFTH  YEAR. 

(a)  May  (c)  July 

(b)  June  (d)  August 

Solid  contours  indicate  deepening  of  the  interface 
and  antir”  Ionic  circulation;  dashed  contours  indicate 
shallowi.  f  cyclonic  circulation . 
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INTERFACE  DEVIATIONS 


FIG.  IS  SOLUTIONS  OF  COMBINED  MIND  AND  INFLOW  FORCING  EXPERIMENT 
DURING  THE  FIFTH  YEAR. 

(a)  September  (c)  November 

( b )  October  (d)  December 

Solid  contours  indicate  deepening  of  the  interface 
and  anticyclonic  circulation;  dashed  contours  indicate 
shallowing  and  cyclonic  circulation. 
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